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Singularities and asymptotic behavior of the 
Tolman-Bondi model * 

Alexander Gromov ^ 


Abstract 

The Bondi formula for calculation of the invariant mass in the 
Tolman- Bondi (TB) model is interprated as a transformation rule on 
the set of co-moving coordinates. The general procedure by which 
the three arbitrary functions of the TB model are determined explic¬ 
itly is presented. The properties of the TB model, produced by the 
transformation rule are studied. 

Two applications are studied; for the falling TB flat model the 
equation of motion of two singularities hypersurfaces are obtained; 
for the expanding TB flat model the dependence of size of area with 
friedmann-like solution on initial conditions is studied in the limit 
t —> +00. 

PACS number(s): 95.30.Sf, 98.65.-r, 98.80.-k 
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1 Introduction 

The Tolman-Bondi (TB) model (Tolman, 1934; Bondi, 1947) describes a 
spherical symmetry self consistency dnst motion in a co-moving coordinates. 
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The Einstein’s eqnations are rednced to the definition of the coefficient gi^i of 
the metrical tensor, the second order partial differential eqnation of motion 
for the fnnction i?(r, t) (the analog of Euler coordinate) and the algebraic 
density law. 

The Tolman-Bondi (TB) model depends on three undetermined functions: 
/(r), F{r) and F(r) of a co-moving coordinate r. One of them, /(r), defines 
the gi i component of the metric and two others are obtained in the process 
of integration of the equation of motion. 

The TB model is defined up to a transformation : r —r. But the TB 
model is also a hydrodynamical problem which requires only two undeter¬ 
mined function to start with; the initial density and velocity distribution. 
In the set of publications (Ruban and Chernin, 1969; Bonnor, 1972; Bonnor, 
1974; Silk 1977; Gurevich and Chernin, 1979) the number of undetermined 
functions decreases by the use of the different transformation of co-moving 
coordinates (see table 1). The common property of these transformations is 
the definition of new coordinates through not more that the only one of three 
undetermined functions of the TB model. 

The main idea of this paper is the interpretation of the the invariant 
mass formula, calculated by Bondi (1947), as the transformation rule from 
an arbitrary co-moving coordinates r to the invariant mass coordinate Mj. 
The Mi coordinate is singled out by its physical sense. The coordinate Mi 
will be called Bondi coordinate and the definition of invariant mass in the 
TB model will be called Bondi transformation. 

Unlike other transformations, the Bondi transformation is defined on the 
couple of undetermined functions (/, E), so that the dependence F = F{f) 
is provided. It is the new way to decrease the number of undetermined 
functions and the way is unique for all meanings of /. 

The Bondi transformation produces also: the dependence of initial con¬ 
dition Ro(^) for fhe function R{r,t) on the function / and initial density 
distribution po{r): Rq = i?o(/;Po(r)); the integral equation for the function 
/ with coefficients po{r) and initial velocity profile Ro{r)] the factorisation on 
the set of couples of functions (/, E); the factorisation on the set of couples 
of initial conditions (i?o(r), po(r)). 

The obtained results are applied to study several properties of the fiat 
TB model. For the falling fiat models the equations of the singularity hy¬ 
persurfaces are obtained. It is shown how the gradient of the initial density 
distribution produces the split of one singularity of FRW model to two sin- 
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author 

4> : r —new coordinate 

Ruban and Chernin,1969 
Bonnor, 1972 

Bonnor, 1972 

Bonnor, 1974 

Silk, 1977 

Gurevich and Chernin 1979 

r —M(r) 
r —>• a Roir) 

/-I 1 3\~1/6 

r ^ a [1 + ar ) 
r ^ a F^/^(r) 
r ^ F{r )/(2 F(r)) 
r —F(r)/4 


Table 1: Examples of co-moving coordinates. M is the mass(in the Newto¬ 
nian theory); r is the arbitrary co-moving coordinate; Ro{r) is the analog of 
Euler radial coordinate; a is the constant; E{r) is the full specihc energy. 


gularities of the TB model. For the expanding TB models the dependence of 
a character of friedmannisation on the initial density distribution is studied. 

First the statement about friedmannisation of the expanding TB models 
has been presented in the article Collins and Hawking (1973). They have 
studied the homogeneous cosmological models but they noted that although 
they work ’... does not prove that there is no open set of inhomogeneous ini¬ 
tial data which gives rise to models that approach homogeneity and isotropy 
... it makes it appear very improbable since one would expect inhomogeneities 
to produce anisotropy rather that isotropy’. Bonnor (1974) has presented the 
counterexample and sowed that ’... an expanding Tolman model with P = 1 
and everywhere non-vanishing density necessarily avolves to the homoge¬ 
neous Einstein- de Sitter model, whatever its initial condition.’ This paper 
studies the same problem from the viewpoint of the results obtained on the 
ground of the Bondi transformation. 

The short review of the TB model is presented in the section The 
interpretation of the Bondi formula for calculation the invariant mass as a 
transformation rule and the dehnition of the functions F and F on the base 
of the transformation rpesents in the section In the section ^ one initial 
condition for the equation of motion and the equation for function / are 
obtained; the partial case: FRW model is discussed. The sections || and || 
are dedicated to the application of the obtained results to the falling and 
expanding flat TB models. The article is ended by the summary. 
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2 The TB Model 


The Tolman-Bondi model (TB) is the simplest spherical symmetry inhomo¬ 
geneous model in general relativity (Tolman, 1934; Bondi, 1947). TB model 
describes pressureless dust motion in a co-moving coordinates in the metric 
(c= 1,^ = 1) 

ds^{r,t) = dt^ - - R^{r,t)dn^, ' = (1) 

where R{r,t) is an analog of Euler coordinate, /(r) is the undetermined 
function of the model, obtained by the integration of (1,4) component of the 
Einstein’s equations. Others equations give the equation of motion 

2 R{r, t)R{r, t) -|- R‘^{r, t) + I — P{r) = 0 (2) 


and the density law: 


8vrp(r, t) = 


dF(r) 1 


(3) 


The equation of motion is solved together with initial conditions: 


^(Bf)L=o = Ro{r), 


R{r,t) 


t=o 


= i?n(r) 


(4) 

(5) 


here Ro{r) and Ro{r) are specified functions. The first integral of the equation 
of motion is 




+ 


F{r) 


( 6 ) 


4i?(r,t)’ 

where F{r) is the second undetermined function. Bondi got the interpre¬ 
tation of two functions / and F on the ground of the comparison of the 
equation (P) with Newtonian theory: ^ may be interpreted as an ana¬ 
log of a full specihc energy Eo{r) and F{r)/4 may be interpreted as an analog 
of an effective gravitating mass. 

The integral of the equation (0) has the form 


±t + F(r) = 


dR 


\JPir) - 1 + 


F{r) ’ 


(7) 


2R 
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where F(r) is the third undetermined function, the sign corresponds to 
an expanding solution and the sign 'corresponds to the falling one. 

The functions /(r), F{r) and F(r) are three undetermined functions of 
the TB model. 


3 The Definition of the Functions F and F 


Ruban and Chernin (1969) have studied the Newtonian analog of the TB 
model in the coordinatisation, dehned by the equality 

r = M. (8) 


It was shown that the use of this equality decreases a number of undetermined 
functions by one unit. 

Bondi showed that in the Newtonian theory invariant mass and gravitat¬ 
ing mass coincide, but in the TB model it is not true (Bondi, 1947). The 
using of the effective gravitating mass coordinatisation given by F/4 has been 
represented in the set of articles (for example, Gurevich and Chernin (1979), 
Silk (1977) and it is also the way to decrease the number of indeterminated 
functions from 3 to 2. 

In this paper the invariant mass Mj, calculated by Bondi (1947): 

r 

Mi(r) = 47r J p{r,t) t) dr = 

0 


47r 



t) dR{r, t) 
/(r) dr 


dr 



F'{r) 

fir) 


dr. 


F'(r) 


dF(r) 

dr 


(9) 


is used as independent co-moving variable. We will name Mj Bondi coor¬ 
dinatisation. The fact that Mj is independent from time means that the 
particle layers do not intersect. The dehnition of invariant mass has a form 
of a transformation rule from r to Mj. This transformation will be called 
Bondi transformation. The identical Bondi transformation 


r = Mj 


( 10 ) 
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produces the relation between two functions / and F\ 

Mi 

F(M,) =4 I /(M,)dM,. (11) 

0 

/^From this point in this paper we will use the Bondi coordinatisation. /^From 
(H) and (1^) it follows that the coordinatisation Mi is singled out by the fact 
that the invariant density is constant: 

p{Mi,t) ^-g{Mi,t) = ( 12 ) 

’ 4 TT 

Now let us study the case of two coordinatisations: an undetermined 
coordinatisation r and M^. Suppose the functions /(r) and F(r) are specihed. 
The Bondi transformation is the example of the continuous transformation 
up to which the coordinatisation is dehned. So, we have three functions /(r), 
F{r) and Mi{r), any two of which dehne the third function. 

Tolman (1934) also uses the third undetermined function F to solve the 
equation (^. The function F(Mj) in the coordinatisation Mj has the form 

F{Mi) = J I ~ ( 13 ) 

y/2(M,)-l + J //(M,)dM, 

We can see now that two functions: F and F depends on the function /. 


4 The Initial Conditions for the TB Model 


Let us suppose that an initial prohle of density is given in the Bondi coordi¬ 
natisation. From (^ it follows that 


R(Mi, 0 ) — Roi^Mi) 


3 

4 TT 


7 fm 

I p(M.,0) 


1/3 


dM, 


(14) 


/^From (|T^) and (|^) it follows that / > 0 and F > 0. Using ([TT|) and (H) in 
the initial moment of time t = 0, together with the initial condition ([l^) , we 
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obtain the following eqnation for the fnnction f{Mi): 

1/3 


1 

2 


3 \ 


1/3 


471 


/ 


•f{Mt)dMi\ 

I 


Mi 


RliMi) - f (M,) + 1 = / /(M,)dM,. 


The solntion of this equation has the form 


(15) 


(16) 


/^From ([I^ it follows that f{Mi) > 0 for every initial density distribution 
Po{Mi) , and from O it follows that F > 0. This nonequality implies 
RliMi) - p{Mi) + 1 > 0. 

The equation should be used as follows: we start from the initial condi¬ 
tions po{Mi) and Ro{Mi) and calculate the function Then we £x the 

function /(Mj) and (|I^ becomes the definition of the set of initial conditions 
which do not change the function /(Mj). 

Substituting (|TT[) and ([T^) in we see that the density law becomes the 
identity in the initial moment of time t = 0. So, the initial conditions of the 
TB model coincide with the initial conditions of the equation of motion (j^). 

In case of FRW model the initial condition ([T^ becomes 


RpRwiMi, 0 ) — 


PFRW {^) 


Mi 1 1/3 

[ /(M,) dm 


In case of / = 1 and p{Mi, 0) = po = const (^) is reduced to 


3H^ = 8 npo, 


where 


H = 


Ro{Mi) 


ISnpo 


FRW 


RoiMi) 

On the other hand, from the assumption that 

Ro{Mi) = H Ro{M,), 

where H = const, from ([T^) it follows that 

Po{Mi, 0) = const. 


= const. 


( 17 ) 

(18) 

(19) 

( 20 ) 
( 21 ) 
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5 Application for Falling Flat Solution 


In this section we will study several properties of the falling flat {p = 1) TB 
model. The solution of the equation of the TB model with A = 0 has been 
obtained by Bonnor (1972, 1974). The flat Bonnor’s solution is obtained by 
the substitution of /^ = 1 into (|]) and the following equation is obtained: 




(dR{M„t) 
\ dt 


2 


2M,. 


( 22 ) 


This equation is studied together with the initial conditions (|T^). The solu¬ 
tion is as follows: 





(23) 


Note here, that from (^) and (0 it follows that the possibility of repre¬ 
sentation of the R{Mi,t) as the multiplication of the two functions, one of 
which depends only on Mj and the second one, which depends only on t, is 
truly only for FRW model: 


RFRwi^i^ 






(24) 


The density law of the flat TB model is obtained by the substitution of 
(^31) into d^) and has the form: 


4:71 p{Mi,t) = J 


p3/2 

Kq 


m - ^ 


/M 


R^{M,) R',{M,)- 


/2Mi 


(25) 


The density has the singularity when one of the parts of the equation in 
parenthesis is equal to zero. Two singlarities are in accordance with two 
difference reasons: the falling of the dust to the cenrte i? = 0 is described by 
the equation 


72 
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RT{M,) 


M, = 0; 


(26) 















the breakdown of the condition of particle layers nonintersecting eqnivalents 
to the eqnation 


AW) RiiMi) - = 0. (27) 

ti(Mj) and t 2 {Mi) are the hypersnrfaces where p{Mi,t) = oo. 

Substituting ([I^ into ( p6D and ( 0 ) we obtain the equations of the hy¬ 
persurfaces: 





Mi 


dMi 

PoiMiY 


(28) 


t2{MY 


po{Mi) 


\ 


Mi 


Mi 

6vr / 
0 


AMi 

Po(Mi) 


/^From ( ]28| ) and (^) it follows 


tYMi)t2{Mi) 


1 

6 7rpo{Mi)' 


In the limit of p{Mi,t) = p{t) we obtain 


tYMi)=t2{Mi) 


1 

V^6irp(«) 


(29) 


(30) 


(31) 


This means that a gradient in the initial prohle of the density splits one 
singularity of the FRW model into two singularities in the TB model. 


6 Application for Expanding Flat Solution 

In this section we apply the results of sections and to the expanding 
TB flat model. 

Bonnor (1974) showed that the set of flat TB models evolve asymptoti¬ 
cally towards a FRW model. This is meant to be a counterexample to the 
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statement of Collins and Hawking (1973) that measure zero of all models can 
evolve towards FRW. 

Using the initial conditions in the form (0) and Bondi coordinatisation 
Mi we will study here the influence of different initial conditions on the 
character of friedmannisation of the Tolman-Bondi flat models. 


6.1 Bonnor’s initial condition 


First we will analyze the initial conditions given by Bonnor (1974). Bonnor 
uses the flat solution of the TB model in the form: 

R(M„f) = 1 [^F{Mi)]^/^[t + f3{Mi)f/\ (32) 


Using (0) together with p 
that 


1 and comparing (p^ and (p^ we hnd out 


R(M„t) 


91/3 




(33) 


where 


Using the initial conditions (0P we obtain 


(34) 


Mi 

1 f dMi 


\\GnMi I p,{Mi) 

Bonnor (1974) studies one particular case of initial conditions 

/3o = /3(0) = 0. 

It follows from (^) that 

/3(0) = ^ 


(36) 


6irpi)(0) 


(36) 


(37) 
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means 


Po(0) = oo. (38) 

At the same time it follows from and ( |36D that 

i?(M„0) = 0, (39) 

and there are no particles in the area Mj >0. It means that the Bonnor’s 
initial condition (^) produces the TB model with delta-like distribution of 
dust: 


Po(M,) =5(M,). (40) 

Bonnor has used the equation (|36|) to prove the statement: an ex¬ 

panding Tolman model with /^ = 1 and everywhere non-vanishing density 
necessarily avolves to the homogeneous Einstein- de Sitter model, whatever 
its initial condition.’ Using the Bonnor transformation, we saw that this 
statement is true only for delta-like initial density distribution. 


6.2 Non delta-like initial conditions 


We will study now how the non delta-like initial condition for density 


p(Mi,0) = pQ^Mi) 


(41) 


change the result. This initial condition produces new asymptotic properties 
of the TB model. 

We will study the following class of initial conditions of the density dis¬ 
tribution: the density monotonically decreases with increase of Mj and 


lim PoiMi) = 0 


(42) 


together with 


Mi 


lim — 

R(Mi)^+oo Mi 


dMj 

po{Mi) 


= +CX). 


The left side of the equation is the average of 1/po(Afj). 
initial conditions 

lim ti(Mj) = -|-cx) 

R(Mi)^+oo 


(43) 

For this class of 

(44) 
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and 


t2{Mi) > ti(Mi), 


(45) 


SO 

lim max{t 2 {Mi),ti{Mi)\ = + 00 . 

R(Mi)^+oo 

The inversion of the eqnations (^) and (pQ]) has the form 
Ml = Mi{t) and M 2 = M 2 {t), 


(46) 


(47) 


and describes two snrfaces in the Tolman space-time with the following prop¬ 
erties: the snrface Mi (t) separates the area 0 < Mi < Mi (f) where the initial 
conditions are forgotten for the fnnction 

R{Mi,t > fi(Mj)) and the area Mi{t) < Mi, where it is not possible to 
neglect the initial conditions; similarly, in the area 0 < Mj < M 2 {t) the ini¬ 
tial conditions are forgotten for R'{Mi,t) and for time t > t 2 {Mi). For the 
velocity dMi/dt we obtain: 


dMi(f) 
hm -^- 

R{Mi)^+oo dt 


0 . 


(48) 


The conditions and (|^) show that for the chosen class of the initial 
conditions both snrfaces Mi(f) and M 2 {t) never arrive in the point where 

Po{Mi) = 0 . 

Let ns calcnlate now the density in the two areas: 

0 < Mcore <'min{Mi{t), M2{t)} (49) 

and 

max{Mi{t), M 2 {t)} < M^heii < + 00 . (50) 

The density law for the expanding model is: 

^ Ovr MMi)+t] [t2iMi)+ty 

where fi(Mj) and t 2 {Mi) are represented by (|^) and (|^) respectiveles. At 
the limit t ^ +00 from (0) it follows that 

pcore{Mi,t) ^ (52) 

ottR 
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and 


Pshell^^i') 



(53) 


The both initial conditions for R{Mi, t) and for R'{Mi, t) are forgotten simul¬ 
taneously in the core in the limit of 


t 3> max{ti{Mi), 


(54) 


Because the both surfaces Mi{t) and M 2 {t) never arrive in to the point 
with pQ{Mi) = 0 the area of friedmannisation of the TB flat expanding model 
is also limited by the nonequality 


Mcore -C min{Mi(t), M 2 (t)}. 


(55) 


7 Summary 

The use of the dehnition of the invariant mass allows to dehne the transforma¬ 
tion rule on the set of co-moving coordinates. The identical transformation 
is studied in this paper. 

The indentical transformation produces the initial conditions for the func¬ 
tion is given by (JT^ and the dehnition of the function /(Mj) is given 

by ([I^). The equation has the solution in the form / = f (pQ^Mi), Ro^Mi)), 
so the function / dehnes the TB model. For hxed TB model, i.e. for hxed 
function /, the equation (|^ becomes the dehnition of the set of equivalent 
couples of initial conditions {po{Mi), Rq^Mi)). 

The hat models are studied together with initial conditions ([I^) . For the 
falling hat models the equations of the singularity hypersurfaces are studied. 
It is shown how the gradient of the initial density distribution produces the 
split of one singularity of FRW model to two singularities of the TB model. 
For the expanding hat model it shown that only delta-like initial density 
distribution produces the FRW model in the limit of t -fcxo. For the 
chosen class of initial conditions non delta-like initial density distribution go 
to the FRW-like core and TB-like shell in that limit. 

The obtained solutions for the falling and expanding hat models are 
stronly depends on the initial conditions but independs from the value of 
function /. It means that every expanding (i.e. non hats) TB models 
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also have the friedmann-like core and nonfriedmann-like shell in the limit 
t +CX). The difference is in the structure of the TB-like shell. The mo¬ 
tion of the singularitiy hypersurfaces (for falling TB models) has a different 
character for different initial density prohles and can brake the toplogy of 
the area of dehnition of the TB model. 
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